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Summary 

A quantum computer has a clear advantage over a classical computer for exhaustive search. The quantum 
mechanical algorithm for exhaustive search was originally derived by using subtle properties of a particular quantum 
mechanical operation called the Walsh-Hadamard (W-H) transform. This paper shows that this algorithm can be 
implemented by replacing the W-H transform by almost any quantum mechanical operation. This leads to several 
new applications where it improves the number of steps by a square-root. It also broadens the scope for implementa- 
tion since it demonstrates quantum mechanical algorithms that can readily adapt to available technology. 

0. Introduction Quantum mechanical systems can be in a superposition of computational states and hence 

simultaneously carry out multiple computations in the same computer. In the last few years there has been extensive 

research on how to use this quantum parallelism to carry out meaningful computations. In any quantum mechanical 

computation the system is initialized to a state that is easy to prepare and caused to evolve unitarily, the answer to the 

computational problem is deduced by a final measurement that projects the system onto a unique state. The ampUtude 

(and hence probabihty) of reaching a specified final state depends on the interference of all paths that take it from the 

initial to the final state - the system is thus very sensitive to any magnitude or phase disturbances that affect any of the 

paths leading to the desired final state. As a result quantum mechanical algorithms are very dehcate and it is believed 

an actual implementation would need elaborate procedures for correcting errors [Err]. 

This paper shows that the quantum search algorithm is surprisingly robust to certain kinds of perturbations. It 
was originally derived by using the W-H transform and appeared to be a consequence of the special properties of this 
transform, this paper shows that similar results are obtained by substituting any unitary transformation in place of the 
W-H transform. Since all quantum mechanical operations are unitary, this means that any quantum mechanical sys- 
tem can be used - all that is needed is a valid quantum mechanical operation and a way of selectively inverting the 
phase of states. Meaningful computation can hence be carried out on the basis of universal properties of quantum 
mechanical operations, this is somewhat similar in spirit to [Neural] where circuit behavior of a certain class of neural 
networks was independent of the precise nature of the nonlinearity in each neuron. 

1. Quantum operations in a quantum computer, the logic circuitry and time steps are essentially classical, 
only the memory bits that hold the variables are in quantum superpositions - these are called qubits. There is a set of 
distinguished computational states in which aU the bits are definite Os or Is. In a quantum mechanical algorithm, the 
quantum computer consisting of a number of qubits, is prepared in some simple initial state, and caused to evolve uni- 
tarily for some time, and then is measured. The algorithm is the design of the unitary evolution of the system. Opera- 
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tions that can be carried out in a controlled way are unitary operations that act on a small number of qubits in each 
step. Two elementary unitary operations presented in this section are the W-H transformation and the selective inver- 
sion of the amplitudes of certain states. 

A basic operation in quantum computing is the operation M performed on a single qubit - this is represented by 



the following matrix: M = 

J2 



1 1 
1 -1 



- the state is transformed into a superposition: | | . Similarly state 1 



is transformed into the superposition \-^, — ^ | . In a system in which the states are described by n qubits (it has 
N = 2 possible states) we can perform the transformation M on each qubit independently in sequence thus chang- 

fl fl 

ing the state of the system. The state transition matrix representing this operation will be of dimension 2x2. Con- 
sider a case when the starting state is one of the 2" basis states, i.e. a state described by a general string of n binary 
digits composed of some Os and some Is. The result of performing the transformation M on each qubit will be a 

n 

superposition of states consisting of all possible n bit binary strings with amplitude of each state being ±2 . This 
transformation is referred to as the W-H transformation [DJ]. This operation (or a closely related operation called the 
Fourier Transformation [Factor]) is one of the things that makes quantum mechanical algorithms more powerful than 
classical algorithms and forms the basis for most significant quantum mechanical algorithms. 

The other transformation that we wiU need is the selective inversion of the phase of the amplitudes in certain 
states. Unlike the W-H transformation and other state transition matrices, the probabihty in each state stays the same 
since the square of the absolute value of the amplitude in each state stays the same. Its realization is particularly 
straightforward; based on [BBHT], we give such a realization. 

Assume that there is a binary function /(x) that is either or 1 . Given a superposition over states x , it is possi- 
ble to design a quantum circuit that will selectively invert the amplitudes in all states where f{x) = 1 . This is 
achieved by appending an anciUa bit, b and considering the quantum circuit that transforms a state \x, b) into 
|x, f(x)XOR d) (such a circuit exists since, as proved in [Reversible], it is possible to design a quantum mechanical 
circuit to evaluate any function/(jc) that can be evaluated classically). If the bit b is initially placed in a superposi- 
tion ^(|0) - 11)) , this circuit will invert the amplitudes precisely in the states for which f{x) = 1 , while leaving 
amplitudes in other states unchanged. 



2. Amplitude amplification a function f{x),x = l , 2 . . . , is given which is known to be non-zero at a 
single value of x , say at x = T - the goal is to find t . If there was no other information about f(x) and one were 
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A'' 

using a classical computer, it is easy to see that on the average it would take — function evaluations to solve this 

problem successfully. However, quantum mechanical systems can explore multiple states simultaneously and there is 
no clear lower bound on how fast this could be done. [BBBV] showed by using subtle arguments about unitary trans- 
forms that it could not be done in fewer than Q.(Jn) steps - subsequently [Search] found an algorithm that took 
exactly 0(JN) steps^ 

The basic idea of [Search] was to consider an A'^ state quantum mechanical system and map each value of x to 
a basis state of the system. The system was initialized so that there was an equal amplitude in each basis state, then by 
a series of unitary operations, the amplitude in the state corresponding to x = T is increased (this state is denoted as 
the |x) state). A measurement is then made due to which the system collapses to a basis state, the observed basis state 
then gives the solution to the problem with a high probability. This algorithm was based on subtle properties of the W- 
H transform. The analysis in this section shows that very similar results are obtained by replacing the W-H transform 
by any arbitrary unitary operation. Some of the consequences of this are presented in the next section. 

Assume that we have at our disposal a unitary operation U and we start with the system in a basis state that is 
easy to prepare, say |y) . If we apply U to |y) , the amplitude of reaching state x is U.^ , if we were to observe the 
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system at this point, the probability of getting the right state would be U \ - according to the notation, y denotes 

I Ml 



the initial state and |x) the target state. It would therefore take Q. 



' 1 ^ 



repetitions of this experiment before a sin- 



gle success. This section shows how it is possible to reach state x in only 0[ , ^ . I steps. This leads to a sizable 



improvement in the number of steps if I U-J « 1 . 

I n 



Denote the unitary operation that inverts the amplitude in a single basis state x by . In matrix notation this is 

the diagonal matrix with all diagonal terms equal to 1 , except the xx term which is -1 - a quantum mechanical 
implementation of this was presented at the end of section 1. 

Consider the following unitary operator: Q = -I U ^I^U - note that since U is unitary, U ^ is equal to the 

adjoint (the complex conjugate of the transpose) of U . We first show that Q preserves the two dimensional vector 

space spanned by: |y) and U V) (note that in the situation of interest, when U is small, these two vectors are 
almost orthogonal). 



1. 0(f(x)) means asymptotically less than a constant times f(x) , 
^(f(x)) means asymptotically greater thsai a constant times f(x) . 
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First consider Qly) . By the definition of Q , this is: -I.,U I^U\y) . Note that |x)(x| , where x is a basis state, 



is an X square matrix all of whose terms are zero, except the xx term which is 1 . Therefore 1^ = 1 - 2|x)(x| & 
7^ = 7- 2|y)<yI, it follows: 

(1) Qlj) = -iI-2\y}{y\)U~\l-2\x}{x\)U\y} = -{I -2\y}{y\)U~^U\y} + 2(1 -2\y}{y\)U~\x}{x\U\y} 
Using the facts: U = I and (yIy) = 1 , it follows that: 



(2) eiy) = ly) + 2(7 - 2|y)<y|)f/ \\x){x\}U\y) . 

-1 * 

SimpUfying the second term of (2) by the following identities: (t| U\y) = U & (yl U |x) = U (as mentioned pre- 

%j zy 

viously, U is unitary and so J7 ^ is equal to the complex conjugate of its transpose): 

(3) eiy) = IY)(1 -4|f/^/) +2f/^^(f/"^|T)) 



Next consider the action of the operator Q on the vector U |x) . Using the definition of Q (i.e. 
Q = -I^U ^ I^U) and carrying out the algebra as in the computation of Q|y) above, this yields: 

(4) e(C/~^|T)) = -7^f/~^7^C/(f/"^|x)) = -7^f/"^7^|x) = I^U~\x) . 
Writing as 1^ = 1- 2|y)<y| and as in (3), <y| f/" V> = f/*^ : 

(5) QiU'^\x)) = t/'^|x) + |yXy|(t/"^|x» = t/"^|x) - 2t/*^|y) . 



1, 



It follows that the operator Q transforms any superposition of the vectors |y) & U |x) into another superpo 



sition of the same two vectors, thus preserving the two dimensional vector space spanned by \y} & U ^ |x) . (3) & (5) 
may be written as: 



(6) Q 



IY> 



(1 -4|f/J ) 2U 



xyl 



xy 



-2U^ 



xy 



lY) 



This yields a process with a periodicity of ^i^j — j as in [BBHT]. If we start with |y), then after -^r^ — , repetitions of 



xyl 



xyl 



Q we get the superposition defined by i7 ^ |x) . From this, with a single application of i7 , we can get |x) . Therefore 
1 



in O 



U 



xyl 



steps, we can start with |y) and reach the target state |x) with certainty. 



The above derivation easily extends to the case when the amplitudes in states, y cfe x , instead of being inverted 
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by & I-, are rotated by arbitrary phases. However, the number of operations required to reach t will be greater. 
Given a choice, it would be clearly better to use the inversion rather than a different phase rotation. Also the analysis 
can be extended to include the case where is replaced by V ^I^V , V is an arbitrary unitary matrix. The analysis 
is the same as before but instead of the operation U , we wiU now have the operation VU . 



3. Examples of quantum mechanical algorithms The interesting feature of this paper is that u 
can be any unitary transformation. Clearly, it can be used to design algorithms where f/ is a transformation in a quan- 

tum computer - this paper give a few such apphcations. The N = 2 states to be searched are represented by n 
qubits. According to the framework of section 2, a search of the A'^ states can be carried out quantum mechanically, if 
we have a unitary operation U which has a finite amplitude i7_„ to go from the starting state Y to the target state T . 

Such a search will take 0\ .} A steps. The following analyses calculate U.^ and thus the number of steps required 

for various algorithms - (i) & (ii) use the W-H transform as U ; (iii) uses a different unitary operation. 

(i) Exhaustive search starting from the state In case the starting state y be the state and the unitary 

transformation U is chosen to be W (the W-H transformation as discussed in section 1), then i7 for any target 

State t is -p . Section 2 gives an algorithm requiring o[ . ^ . \ i.e. OiJN) steps. This algorithm would carry out 

repeated operations of Q, with U ^ = U = W , Q becomes Q = -I-WI-W; the operation sequence is hence: 

^ 

. . .{-I -W I -W){-I -W I -W)i-I -W I -W) . . . By rearranging parentheses and shifting minus signs, this may be seen to 

X ^ ^ 

consist of alternating repetitions of -WI-W & . 

The operation sequence -WI-W is simply the inversion about average operation [Search]. To see this, write /- 

as / - 2|0)<0| . Therefore for any superposition \x) : -WI^W\x) = -W(I - 2|0)<0|)W|jc) = - \x) + 2W\b){b\W\x) . It 

is easily seen that W|0)(0| W\x) is another vector each of whose components is the same and equal to A where 

N 

A = — ^ X. (the average value of all components). Therefore the i''^ component of -WI-W\x) is simply: 
/ = 1 

(-X. + 2 A) . This may be written as A + (A - x.) , i.e. each component is as much above (below) the average as it was 
initially below (above) the average, which is the inversion about average [Search]. 
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(ii) Exhaustive search starting from an arbitrary state For the W-H transform, U^^ between any pair 
of states y & X is ±-p . We can start with any state y and the analysis of section 2, gives us an algorithm to reach x 

in oi .^ , ], i.e. 0{4n) , iterations. Therefore instead of starting with the state, as in (i), we could equally well 

I 'cyl 

start with any state y, and repeatedly apply the operation sequence Q = -I^WI W to obtain an equally efficient 

y X 

0(JN) algorithm. The dynamics is similar to (i); however, there is no longer the convenient inversion about average 
interpretation. 

(iii) Search when an item near the desired state is known 

Problem: Assume that an n bit word is specified - the desired word differs from this in exactly k bits. 
Solution The effect of this constraint is to reduce the size of the solution space. One way of making use of this con- 
straint, would be to map this to another problem which exhaustively searched the reduced space using (i) or (ii). How- 
ever, such a mapping would involve additional overhead. This section presents a different approach which also carries 
over to more complicated situations as discussed in [Qntappl]. 

Instead of choosing U as the W-H transform, as in (i) & (ii), in this algorithm U is tailored to the problem 
under consideration. The starting state y is chosen to be the specified word. The operation U consists of the transfor- 



mation 



4 - 1 ^ 



, appUed to each of the n qubits (a is a variable parameter yet to be determined) - note that 



n- k k 



n 2 fi^2 



if a is 2 , we obtain the W-H transform of section 1. Calculating U^^ it follows that \U^j^ ~ \ a J Voc 

flf fX Tt k k ft k 

this is maximized when a is chosen to be 7 ; then log It/ I = -log -log— — . The analysis of section 2 can 

k I mI 2 n 2 K 

now be used - as in (i) & (ii), this consists of repeating the sequence of operations -^^U ^I^U , tiinss. 



Tyi 



The size of the space being searched in this problem is , which is equal to ^—r- . Using Stirhng's 

\kJ n-klk] 



approximation: logn! = nlogn - n , from this it follows that log , = n\og - k\og , comparing this to the 

\kj n-k n-k 

number of steps required by the algorithm, we find that the number of steps in this algorithm, as in (i) & (ii), varies 
approximately as the square-root of the size of the solution space being searched. 
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4. General quantum mechanical algorithms The framework described in this paper can be used 

to enhance the results of any quantum mechanical algorithm. Assume there is a quantum mechanical algorithm Q due 

to which there is a finite amplitude 2_„ for transitions from the starting state y to the target state x . The probability 

'•I 



of being in the state x is hence \Q^A - it will therefore take Q. 

I Ml 



^ 1 ^ 



repetitions of Q to get a single observation 



of state X. Since the quantum mechanical algorithm 2 is a sequence of r] elementary unitary operations: 

Q^Q2---Qy^^^^^^^^^^^^^^^^^^^ ^^^^^^oxm^^^on. Mso, Q ^ = Q^Q^_y.Q^ Q ^ is given by a sequence 

of the adjoints of the elementary unitary operations in the opposite order and can hence be synthesized. Applying the 
framework of section 2, it follows that by starting with the system in the s state and repeating the sequence of opera- 
tions: -I^Q ^/ 2 , 0[ .^ ? \ times followed by a single application of Q , it is possible to obtain the state x with 



y x^' \\Q , 



certainty. 



5. Sensitivity in order to achieve isolation, quantum mechanical computers generally have to be designed to be 
microscopic - however it is extremely difficult to exert precise control over naicroscopic individual entities. As a 
result, a serious problem in implementing quantum mechanical computers is their extreme sensitivity to perturba- 
tions. This paper synthesizes algorithms in terms of unitary matrices - as shown in sections 3 & 4, this framework can 
always be specialized to a quantum computer based on qubits; however, it can also be applied directly to more physi- 
cal situations, hopefully reducing the need for error correction [Err] . 

For example, consider a hypothetical implementation of the quantum search algorithm where the qubits are the 
spin states of electrons and the W-H transform is achieved by a pulsed external magnetic field. The results of section 
2 & 3 tell us that it does not significantly alter the working of the algorithm if the axes of the magnets, or the periods 
of the pulse are slightly perturbed from what is required for the W-H transform. Any unitary transform, U , close to 

the W-H transform, will work provided both U & U ^ are consistently appUed as specified. 

6. Limitation As shown in [BBHT], it is possible to express several important computer science problems in 
such a way so that a quantum computer could solve them efficiently by an exhaustive search. Even in physics, several 
important problems can be looked upon as searches of domains. Many spectroscopic analyses are essentially searches 
- a rather dramatic example of a recent search was that for the top quark. The framework of this paper could equally 
well be used here. All that is needed is a means to repeatedly apply a specified Hamiltonian that produces various 

12 

phase inversions and state transitions. For example, it took about 10 repetitions of a certain experiment, consisting 
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of interacting a proton and antiproton at high energies, to obtain 12 observations of the top quark [Quark]. Denoting 

2 

the desired state with the top quark by |x) and the initial proton-antiproton state by ly) , it implies that t/.,, is 

I ^yi 

-12 -1 
approximately 12 x 10 . Therefore if it were possible to apply the operation -I^U IM repetitively m times, it 

2 2 
would boost the success probability by approximately m (assuming m to be a small number), and it would take m 

fewer experiments; in case it were possible to apply the operations -I^U ^I^U, about a million times, one could 
achieve success in a single experiment! 

In principle it is possible to synthesize U ^ for any unitary operation, U , since the adjoint of a unitary opera- 
tor is unitary and can hence be synthesized quantum mechanically. For controlled operations on qubits, synthesizing 
the adjoint is no harder than synthesizing the original operation as discussed in section 4. However, the adjoint of the 
time evolution operation is the reversed evolution operation - this may not be easy to synthesize when the states are 
non-degenerate and there is significant time evolution. This is especially true if the time-evolution is due to the inter- 
nal dynamics of the system. That is the main reason this procedure, at least in its present form, could not be used to 
isolate the top quark! 

This paper shows that the W-H transform of search based algorithms can be replaced by any unitary operation, 
provided the selective inversion operation is carried out precisely. As mentioned at the end of section 2, even the 
selective inversion can be replaced by a selective phase shift, provided it only affected the concemed states. Also, as 

mentioned in section 2, the analysis stays virtually the same if is replaced by V ^IxV with V unitary. 

Another limitation is that the framework of section 2 demands that U & U ^ stay the same at all time steps. 
What happens if there are small perturbations in these? It seems plausible that these will not create much impact if 
they are small and average out to zero; however, that is something still to be proved. 

7. Conclusion Designing a useful quantum computer has been a daunting task for two reasons. First, because 

the physics to implement this is different from what most known devices use and so it is not clear what its structure 

should be like. The second reason is that once such a computer is built, few appUcations for this are known where it 

will have a clear advantage over existing computers. This paper has given a general framework for the synthesis of a 

category of algorithms where the quantum computer would have an advantage. It is expected that this formalism will 

also be useful in the physical design of quantum computers, since it demonstrates that quantum algorithms can be 

implemented through general properties of unitary transformations and can thus adapt to available technology. 

The advantage of the bit representation is that the number of states that can be represented is exponential in the 
number of bits. After the success of classical digital computers, most quantum mechanical schemes, based them- 
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selves on qubits. It was shown that it was possible, in principle, to synthesize gates that operate on qubits [Bit]. How- 
ever independent qubits are not such a natural entity in the real quantum world as bits are in the classical world. They 
are very delicate entities and the slightest perturbation affects the entire computation. The distinctive feature of this 
paper is that it is based on general properties of unitary transforms and can hence be applied to a variety of situations 
including, but not limited to, qubit operations. 

8. Acknowledgments Thanks to Norm Margolus & CharUe Bennett for spending the time and effort to com- 
ment on several versions of this paper. 



9. References 

[BBHT] M. Boyer, G. Brassard, P. Hoyer & A. Tapp, Tight bounds on quantum searching, Proc, PhysComp 

1996 (lanl e-print quant-ph/9605034). 
[BBBV] C. H. Bennett, E. Bernstein, G. Brassard & U.Vazirani, Strengths and weaknesses of quantum com- 
puting, SIAM Journal on Computing, pp. 1510-1524, vol. 26, no. 5, Oct 1997. 
[Bit] S. Lloyd, A potentially realizable quantum computer. Science, vol. 261, pp. 1569-1571, 1993. 

[DJ] D. Deutsch & R. Josza, Rapid solution of problems by quantum computation, Proc. Royal Society of 

London, A400, 1992, pp. 73-90. 
[Err] P. W. Shor, Fault-tolerant quantum computation, Proc. 35th Annual Symposium on Fundamentals of 

Computer Science (FOCS), 1996. 
[Factor] P. Shor, Algorithms for quantum computation: discrete logarithms and factoring. Proceedings, 35th 

Annual Symposium on Fundamentals of Computer Science (FOCS), 1994, pp. 124-134. 
[Neural] J. Hopfield, Neural networks & physical systems with emergent collective computational abilities 

Proceedings National Academy of Sciences, vol. 79, p. 2554, 1982. 
[Qntappl] L. Grover, A framework for fast quantum mechanical algorithms, lanl e-print quant-ph/97 1 1043 . 
[Quark] T. M. Liss & P. L. Tipton, The discovery of the top quark. Scientific American, Sep. '97, pp. 54-59. 

[Reversible] C. H. Bennett, Space-Time trade-offs in reversible computing, SIAM Journal of Computing, vol. 18, 

pages 766-776, 1989. 

[Search] L. K. Grover, Quantum Mechanics helps in searching for a needle in a haystack, Phys. Rev. Letters, 
vol. 78(2), 1997, pages 325-328 (lanl e-print quant-ph/9605043). 



9 



